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This Chapter deals with single layers of carbon (graphene) and hexagonal 
boron nitride on transition metal surfaces. The transition metal substrates 
take the role of the support and allow due to their catalytic activity the 
growth of perfect layers by means of chemical vapor deposition. The layers 
are sp^ hybridized honeycomb networks with strong in plane a and weaker tt 
bonds to the substrate and to the adsorbates. This hierarchy in bond strength 
causes anisotropic elastic properties, where the sp^ layers are stiff in plane 
and soft out of plane. A corrugation of these layers imposes a third hierarchy 
level in bond energies, with lateral bonding to molecular objects with sizes 
between 1 and 5 nanometer. This extra bond energies are in the range of 
thermal energies fc^T at room temperature and are particularly interesting 
for nanotechnology. The concomitant template function will be discussed 
in Section 2.6.2. The peculiar bond hierarchy also imposes intercalation as 
another property of sp^ layer systems (Section 2.5). Last but not least sp^ 
layer systems are particularly robust, i.e. survive immersion into liquids [1], 
which is a promise for sp^ layers being useful outside ultra high vacuum. 

The Chapter shortly recalls the synthesis, describes the atomic and elec- 
tronic structure, is followed by a discussion of properties like intercalation 
and the use of sp^ layers on metals as tunneling junctions or as templates. 
The Sections are divided into subsections along the sketch in Figure 1, i.e. 
for flat and corrugated layers. Since there are flat and corrugated layers for 
graphene (g) as well as hexagonal boron nitride (h-BN) the similarities and 
differences between C-C and B-N are discussed in every subsection. The 
Chapter ends with an Appendix that summarizes the basics of atomic and 
electronic structure of honeycomb lattices. 

Of course the Chapter does not cover all aspects of sp^ single layers. Top- 
ics like free standing layers [2], edge structures of ribbons [3, 4], topological 
defects [5], or mechanical and chemical properties were not covered. 



3 



a) b) 

Molecule 
sp^ layer 

Metal 

Figure 1: Schematic side view of a single layer on top of a metal. The layer 
changes the properties of adsorbed molecules as compared to the pristine 
surface, a) Flat layer, where molecule-molecule interactions are present (see 
Figure 4 and 18. b) Corrugated layer, which is a template on which single 
molecules can be laterally isolated (see Figure 7). 





1 Single layer systems 

A single layer of an adsorbate strongly influences the physical and chemical 
properties of a surface. Sticking and bonding of atoms and molecules may 
change by orders of magnitude, as well as the charge transport properties 
across and parallel to the interface. 

There are many single layer systems like e.g. graphite/graphene [6, 7], 
hexagonal boron nitride [8, 9], boron carbides [10], molybdenum disulflde 
[11], sodium chlorides [12, 13], or aluminium oxide [14], copper nitride [15, 16] 
to name a few. In order to decide whether single layers are "dielectric" or 
"metallic" the electronic structure at the Fermi level has to be studied, where 
a metallic layer introduces new bands at the Fermi energy, while a dielectric 
layer does not. 

Figure 1 shows a schematic view of the single layer systems, with ad- 
sorbed molecules. The beneath conduction electrons of the substrate still 
have tunneling contact to the molecule, though the bonding is much weaker 
compared to the bare metal. If the single layer is flat (Figure 1 a)) ad- 
sorbed molecules may touch, while in corrugated layer systems (Figure 1 b)) 
molecules may be trapped also laterally on the surface. The corrugation (see 
Secton 2.2.2) results from the epitaxial stress due to the lattice misflt be- 
tween the substrate and the sp^ layer. This results in lateral superstructures 
with nanometer periodicities and corrugations in the 0.1 nm range. For de- 
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velopments in nanotechnology it is useful to have single layer systems which 
are inert, remain clean at ambient conditions and are stable up to high tem- 
peratures. In this field sp^ hybridized graphene and hexagonal boron nitride 
are outstanding examples. These two intimately related model systems, are 
discussed in more detail. Both are sp^ hybridized layers with about the same 
lattice constant, but, while on most transition metals graphene is metallic, 
h-BN is an insulator. 
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2 sp^ single layers 



In this Chapter we restrict ourselves to the sp'^ hybridized single layer systems 
of graphene and hexagonal boron nitride. Except for some considerations in 
the Appendix we always deal with layers on transition metal supports. The 
comparison between graphene and boron nitride is particularly helpful for 
the understanding of their functionality. They are siblings with similarities 
and differences. Both form a sp'^ honeycomb network with similar lattice 
constant (0.25 nm), with strong a in plane bonds and soft tt bonds to the 
substrate and the adsorbates. Depending on the substrate, both form flat 
or corrugated overlayers. The corrugation is a consequence of the lattice 
mismatch and the anisotropic bonding of the atoms in the honeycomb with 
the substrate, and leads to the peculiar functionality of molecular trapping 
[17]. Graphene and h-BN differ in their atomic structure: In the case of 
h-BN two different atoms constitute the honeycomb unit cell, while the two 
carbon atoms are equivalent within the honeycomb lattice. This causes most 
graphene overlayers to be metallic, while the h-BN layers are insulators. It 
is furthermore the reason for an " inverted topography" if corrugation occurs 
(see Figure 8). 

The theoretical description of sp'^ single layers [18] is easier than the 
experimental realisation [19]. The flrst experimental reports on sp'^ layer 
production on supports dates back to early times of surface science [20] and 
has sometimes been considered to be an annoyance since they poison catalysts 
[6] . In the field of ultra-thin epitaxial films of graphene and hexagonal boron 
nitride on solid surfaces the work of the Oshima group has to be highlighted 
[21]. They studied the production of sp'^ layer systems systematically and 
carefully characterized their electronic and vibronic structures. Today this 
work is the starting point for numerous ongoing investigations, also because 
of the rise of graphene [2] i.e. the discovery that devices containing free- 
standing single layer devices may be realized. Therefore there is a strong 
demand for production routes for these materials that have, somehow, to 
start on a solid surface. 

2.1 Synthesis 

There are several roads to the production of sp'^ single layers. Here the 
Chemical Vapor Deposition (CVD) processes that base on the reactive ad- 
sorption of precursor molecules from the gas phase onto the substrate, and 
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the segregation method, where the constituents diffuse from the bulk to the 
surface, are briefly summarized. 

The sp^ layers may also be grown on substrates without C^y symmetry 
as e.g. present for the (111) surfaces of face centered cubic fee crystals. 
Examples are the growth of h-BN on Mo(llO) [22] or on Pd(llO) [23]. 

Although this chapter deals with single sp^ layers only, it could be desir- 
able to grow multilayers. There is e.g. one report on the growth of graphene 
on h-BN/Ni(lll) [24], and for the segregation approach it seems to be eas- 
ier to grow multilayers, as it was shown for silicon carbide (SiC) substrates 
[25], since this method has not to overcome the low sticking probability of 
precursor gases on complete sp^ layers. 

2.1.1 Chemical Vapor Deposition (CVD) 

Chemical Vapor Deposition processes comprise the adsorption and reaction 
of the educts i.e. precursor molecules on the surface where a new material 
shall grow. The process often involves the cracking or decomposition of the 
precursor molecules and a partial release of products into the gas phase. 

The first h-BN single layers have been synthesized on Ru(OOOl) and 
Pt(lll) surfaces by CVD of benzene like borazine (HBNH)^ [8]. The pro- 
cess comprises the hydrogen abstraction from the borazine molecules, the 
assembly of hexagonal boron nitride and the desorption of H2 gas. The h- 
BN growth rate drops after the formation of the first layer by several orders 
of magnitude. This has the practical benefit that it is easy to prepare single 
layers. In Figure 2 the growth of h-BN on Ni(lll) is shown as a function of 
the exposure to the precursor molecules. Clearly, the growth rate drops by 
more than two orders of magnitude after the completion of the first layer. 
This behavior is quite general for sp'^ layer systems and similar growth behav- 
ior is expected for graphene layers. The drop in growth rate is presumably 
due to the much lower sticking probability of borazine after the completion 
of the first h-BN layer, and the much lower catalytic activity of h-BN com- 
pared to clean transition metals. However, not much is known on the details 
of the growth process of h-BN. From the study of h-BN island morphologies 
on Ni(lll) it was suggested that the borazine BN six-ring is opened during 
the self assembly process [26, 27]. For h-BN layer formation, also trichlorob- 
orazine (C1BNH)3 [27], and diborane (B2H6) ammonia (NH3) gas mixtures 
[28] were successfully used. 

For the formation of graphene layers many different precursors have been 
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Figure 2: Growth of a h-BN layer on Ni(lll): Coverage as a function of 
borazine exposure. 1 Langmuir (L)= 10~^Torr • s. Note the drop of growth 
rate after the completion of the first layer. By courtesy of W, Auwarter. 



used, as e.g. ethene {C2H4) [29], acetylene {C2H2) [24], or propylene (CsHq) 
[30]. It turns out that almost all hydrocarbons react in non oxidizing en- 
vironments on transition metal surfaces to graphene, though there are also 
growth conditions where diamond films grow [31]. 



2.1.2 Segregation 

An alternative way to the above mentioned CVD processes is the use of 
segregation. Here the educts are diluted in the bulk of the substrate. If 
the substrate is heated, they start to diffuse and eventually meet the surface, 
where they have a larger binding energy than in the bulk. At the surface they 
may then react to the new material. If the substrate temperature exceeds the 
stability limit of the sp'^ layer, back dissolution into the bulk or desorption 
may occur. Well known examples are e.g. the formation of graphene on 
Ni(lll) [6, 32], Ru(OOOl) [33], SiC [25] or the formation of BC3 layers on 
NbB2 [10]. 



2.2 Atomic structure 

The atomic structure of the sp'^ layer systems is fairly well understood. It 
bases on a strong sp'^ hybridized in plane bonded honeycomb lattice and 
a, relative to these a bonds, weak tt bonding to the substrate. The tt (p^) 
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Figure 3: Calculated BN binding energies of h-BN on various transition- 
metals. The h-BN is strained to the lattice constant of the substrate. The 
bond strength correlates with the d-band occupancy, where the bonding is 
strongest for the 4d metals. The dashed lines are guides to the eyes. Results 
within the Local Density Approximation (LDA) are shown. Data from [34]. 



bonding depends on the registry to the substrate atoms, where the layer- 
substrate {pz — dz'2) hybridization causes a tendency for lateral lock-in of the 
overlayer atoms to the substrate atoms. For systems, where substrates have 
the same symmetry (Ca^), as the sp^ honeycombs, the lattice mismatch M 
is defined as: 

^suh 

where aovi is the lattice constant of the overlayer, and a sub that of the 
substrate. Often the sign of the lattice mismatch is not indicated, and the 
absolute value of the difference between the substrate and the overlayer lattice 
constant is used. Then it has to be explicitly said whether the mismatch 
induces compressive (+) or tensile (-) stress in the overlayer, or vice versa, 
tensile (+) or compressive (-) stress in the substrate. The mismatch plays 
a key role in the epitaxy of the sp^ layers. Of course, the bonding to the 
substrate also depends on the substrate type, and in turn may lead to the 
formation of corrugated super structures with a peculiar template function. 

Figure 3 shows the calculated BN bond strength on different transition 
metals (TM) [34]. The calculations were performed for (1x1) unit cells, where 
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the h-BN was strained to the lattice constant of the corresponding substrate, 
with nitrogen in on top position. The BN bond energy is determined from the 
difference of the energy of the (1x1) h-BN/TM system and the strained h-BN 
+ TM system. The trend indicates the importance of the d-band occupancy 
of the substrate and it is e.g. similar to the dissociative adsorption energies 
of ammonia NH^ on transition metals [35]. Since the lock-in energy, i.e. 
the energy that has to be paid when the nitrogen atoms are moved laterally 
away from the on top sites, is expected to scale with the adsorption energy. 
Figure 3 also rationalizes e.g. why h-BN/Pd(lll) is less corrugated than 
h-BN/Rh(lll). 

This section is divided into two parts, where we distinguish flat and corru- 
gated layers. Flat means that the same types of atoms have the same height 
above the substrate, which has only to be expected for (1x1) unit cells like 
that in the h-BN/Ni(lll) system. Corrugation occurs due to lattice mis- 
match and lock-in energy gain, i.e. due to the formation of dislocations. 

2.2.1 Flat layers: Domain boundaries 

The Ni(lll) substrate has Csv symmetry and a very small lattice mismatch of 
+0.4 % to the sp^ layers. On Ni(lll) h-BN forms perfect single layers. Figure 
4 shows scanning tunneling microscopy (STM) images from h-BN/Ni(lll). 
The layers appear to be defect free and flat. The STM may resolve the 
nitrogen and the boron sublattices, where a Tersoff Haman calculation in- 
dicated that the nitrogen atoms map brighter than the boron atoms [36]. 
The production process of the h-BN layers (see Section 2.1) also leads to the 
formation of larger terraces compared to uncovered Ni(lll), where widths 
of 200 nm are easily obtained. The stability of the (111) facet was also ob- 
served in experiments with stepped Ni(755) [37] and Ni(223) surfaces, where 
the miscut of these surfaces relative to the [111] direction lead to large (111) 
facets and step bunches. 

For the case of the h-BN/Ni(lll) system the mismatch M is +0.4%, 
i.e. the h-BN is laterally weakly compressed. This small mismatch leads 
to (1x1) unit cells and atomically flat layers. The atomic structure of h- 
BN/Ni(lll) is well understood and there is good agreement between exper- 
iment [38, 9, 39] and theory [36, 40, 41]. Figure 5 shows six (1x1) configura- 
tions for h-BN within the Ni(lll) unit cell. The nomenclature of the struc- 
ture (B,N)=(fcc,top) indicates that boron is sitting on the fee site, i.e. on a 
site where no atom is found in the second nickel layer, and where the nitrogen 
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Figure 4: Constant current scanning tunneling microscopy (STM) images of 
a flat layer of hexagonal boron nitride on Ni(lll). a) 30x30 nm. The grey 
scales indicate 3 terraces with different heights. Large defect free terraces 
form, b) 3x3 nm. The nitrogen (bright) and the boron (grey) sublattices are 
resolved with different topographical contrast. From [9]. 



atom sits on top of the atom in the flrst nickel layer. Theory found two stable 
structures (B,N)=(fcc,top) and (hep, top) only. In both cases nitrogen is on 
top [36]. The (fee, top) structure has the lowest energy, which was consistent 
with the published experimental structure determinations. The calculated 
energy difference between the structure with boron on fee and on hep dif- 
fers only by 9 meV, which is reasonable since it indicates interaction of the 
overlayer with the second nickel layer. The 9 meV are, however, one order of 
magnitude smaller than the thermal energies ksT during the synthesis, and 
from this viewpoint it is not clear why pure (B,N) = (fcc,top) single domain 
can be grown. Inspection of Figure 5 shows that the (fee, top) structure is a 
translation of the (top, hep) or (hep, fee). A transformation between (fee, top) 
and (hep, top) involves, a rotation, or a permutation of B and N and a transla- 
tion. The growth of the h-BN layers on Ni(lll) proceeds via the nucleation 
and growth of triangular islands that are separated by distances between 
10 and 100 nm [26]. Therefore, if (fee, top) and (hep, top) nucleation seeds 
form, two domains are expected since the change of orientation of islands 
with sizes larger than 10 nm costs too much energy [26]. The concomitant 
domain boundaries (see Figure 16 a)) have interesting functionalities. They 
act as collectors for intercalating atoms, and clusters like to grow on these 
lines. Also, it is expected that such domain boundaries are model systems 
for sp'^ edges. Unfortunately, the (fcc,top):(hcp,top) ratio is not yet under 
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Figure 5: Possible (1x1) configurations for one-monolayer h-BN/Ni(lll). 
The top, fee hollow, and hep hollow sites are eonsidered for the position of 
the boron and nitrogen atoms, respectively. The (B,N) = (fcc,top) registry has 
the strongest bond, while calculations predict 9 meV lower binding energy 
for the (B,N)=(hcp,top) configuration. From [36]. 
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full experimental control. But it is e.g. known that the omission of the oxy- 
gen treatment of the Ni(lll) surface before h-BN growth causes two domain 
systems [26]. 

The vertical position of the nitrogen and the boron atoms is not the same 
within the (1x1) unit cell. This local corrugation or buckling, where nitrogen 
is the outmost atom, and boron sits closer to the first Ni plane, was taken as 
an indication of the compressive stress on the h-BN in the h-BN/Ni(lll) sys- 
tem [42, 9]. The comprehensive density functional theory study of Laskowski 
et al. [34], indicated however that this buckling, i.e. height difference be- 
tween boron and nitrogen persists also for systems with tensile stress in the 
overlayer, where the substrate lattice constant is larger than that of the h- 
BN. B is closer to the first substrate plane for all investigated cases. This is 
a consequence of the bonding to the substrate, where the boron atoms are 
attracted and the nitrogen atoms are repelled from the surface [34]. 

For the corresponding graphene Ni system the same structure i.e the 
(CA,CB)=(fcc,top) which is equivalent to the (top, fee) configuration has been 
singled out against the (fee, hep) structure [38]. There are no reports on 
(top,fcc)/(top,hcp) domain boundaries, as observed for the h-BN/Ni(lll) 
case [26]. For g/Ir(lll), which belongs to the corrugated layer systems, dis- 
location lines that are terminated by heptagon-pentagon defects were found 
[5]. Of course these kinds of defects were less likely for BN, since this would 
imply energetically unfaforable N-N or B-B bonds in the BN network. 

2.2.2 Corrugated layers: Moire and dislocation networks 

When the lattice mismatch M of the laterally rigid sp^ networks exceeds a 
critical value, super structures with large lattice constants are formed. If the 
lattice of the overlayer and the substrate are rigid and parallel, the super 
structure lattice constant gets aovi/\M\^ where aovi is the 1x1 lattice constant 
of graphene or h-BN 0.25 nm). 

For rigid sp^ layers, i.e. if there would be no lateral lock-in energy avail- 
able, we expect flat floating layers, reminiscent to incommensurate moire 
patterns without a directional lock-in of the structures (see Figure 6a)). For 
the case of h-BN/Pd(lll) such a tendency to form moire type patterns, also 
without a preferential lock-in to a substrate direction, were found [44]. These 
h-BN films have the signature of the electronic structure of flat single layers 
(see Section 2.4). The lock-in energy is the energy that an epitaxial system 
gains on top of the average adhesion energy if the overlayer locks into pre- 
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Figure 6: Schematic view of a moire and a dislocation in an overlayer system, 
a) in a moire pattern ^ (n + 1) sp^ units fit on n substrate units, and 
the directions of the substrate and the adsorbate lattice do not necessarily 
coincide, b) in a dislocation network n + 1 sp'^ units coincide with n substrate 
units. Here the lock-in energy is large for on-top and weak or repulsive for 
bridge sites, where the dislocation evolves. 



Figure 7: Relief view of a constant current scanning tunneling microscopy 
image of a corrugated layer boron nitride (nanomesh) (30x30 nm^, 1^=2.5 
nA, V5=-l V). This nanostructure with 3.2±0.1 nm periodicity consists of 
two distinct areas: the wires which are 1.2+0.2 nm broad and the holes with 
a diameter of 2.0±0.2 nm. The corrugation is 0.07±0.02 nm. From [43]. 
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ferred bonding sites. It involves the formation of commensurate coincidence 
lattices between the substrate and the adsorbate layer with dislocations (see 
Figure 6b)). The lock-in energy is expected to be proportional to the bond 
energy shown in Figure 3. For the case of /i-BN/Rh(lll) the lattice mis- 
match (Eqn. 1) is -6.7 %, and 13 x 13 BN units coincide with 12 x 12 Rh 
units [45, 46]. With the room temperature lattice constants of h-BN and Rh 
this leads to a residual compression of the 13 h-BN units by 0.9 %. Figure 7 
shows a relief- view of this super structure coined "h-BN nanomesh" [45]. It 
has a super cell with a lattice constant of 3.2 nm i.e. (12x12) Rh(lll) units 
and displays the peculiar topography of a mesh with "wires" and "holes" 
or "pores". It turned out that this structure is a corrugated single layer 
of hexagonal boron nitride with two electronically distinct regions that are 
related to the topography [47, 48]. This structure also has the ability of 
trapping single molecules in its holes (see Section 2.6.2). The accompanying 
variation of the local coordination of the substrate and the adsorbate atoms 
divides the unit cells into regions with different registries. In reminiscence 
to Figure 5 the notation (B,N)=(top,hcp) refers to the local configuration, 
where a B atom sits on top of the substrate atom in the first layer and N on 
top of the hexagonal close packed (hep) site that is on top of the substrate 
atom in the second layer. Again 3 regions can be distinguished with atoms 
in (fcc,top), (top, hep) and (hep, fee) configurations (see Figure 8). A force 
field theory approach indicated that the (B,N)=(fcc,top) sites correspond to 
the tightly bound regions of the h-BN layers, i.e. the holes, while the weakly 
bound or even repelled regions correspond to the wires [47]. The corrugation, 
i.e. height difference of the h-BN layer from the top of the substrate is, in ac- 
cordance between experiment and theory about 0.05 nm. This is sufficient to 
produce a distinct functionality (see Section 2.6.2) [48]. For h-BN/Ru(0001) 
a structure very similar to that of h-BN/Rh(lll) was found [49]. 

For the graphene case the situation is related, though not identical. The 
difference lies in the fact that the base in the 1x1 unit cell of free-standing 
graphene consists of two identical carbon atoms and C^, while that of 
h-BN does not. and become only distinguishable by the local coor- 
dination to the substrate. In g/Ru the local (fee, top) and (top, hep) coordi- 
nation leads to close contact between the (Ca,Cb) atoms and the substrate 
[51] while (B,N) is strongly interacting only in the (fcc,top) coordination 
[34]. As a result, twice as many atoms are bound in strongly interacting 
regions in g/Ru when compared to /i-BN/Ru. In reminiscence to morpholog- 
ical terms the strongly bound regions of g/Ru(0001) are called valley (V) and 
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Figure 8: Views of the height modulated graphene (g) and /i-BN nanomesh 
(BN) on Ru(OOOl), M and V denote mounds (hep, fee) and vaUeys of the 
graphene, H and W holes (fee, top) and wires of the h-BN nanomesh. The six 
ball model panels illustrate the three different regions ((fee, top), (top, hep) 
and (hcp,fcc)), which can be distinguished in both systems [50]. 
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the weakly bound region with the (0^,0^) atoms on (hep, fee) sites mounds 
(M). The fact that (top, hep) leads to strong bonding for graphene but weak 
bonding for /i-BN gives rise to an inverted topography of the two layers: a 
connected network of strongly bound regions for graphene (valleys) and a 
connected network of weakly bound regions for /i-BN (wires). Also for the 
graphene case "moire"-type superstructures were found, where g/Ir(lll) is 
the best studied so far [7]. 

2.3 Electronic structure I: Work function 
2.3.1 Flat layers: Vertical polarization 

The work function, i.e. the minimum energy required to remove an electron 
from a solid, is material dependent. There are excellent reviews on the topic 
[52, 53]. For our purpose, where we want to discuss the electric fields near 
the surface, it is sufficient to recall the Helmholtz equation that relates the 
work function $ of a ffat surface with vertical electric dipoles: 

$=-A^a-p (2) 
^0 

where e is the elementary charge, eo the permittivity and A^^ the areal 
density of the dipoles p. If the dipoles are assigned to the atoms, we get e.g. 
for Ni(lll) with a work function of 5.2 eV and an in plane lattice constant 
of 0.25 nm a dipole of 0.7 Debye (1 Debye= 3.34 -lO-^^ Cm). It has to 
be said that this is the classical view of the work function and all quantum 
mechanical effects may be incorporated empirically into Equation 2 in using 
an effective dipole. 

Now the inffuence of overlayers shall be analyzed. If a single layer of a 
medium is placed on top of the substrate, the electric fields in the surface 
dipole polarize the layer i.e. decreases the work function by A$5 = -^Na -Pind 
by screening. The induced dipole pind is proportional to the electric field 
perpendicular to the surfaceE"^: pind = <^ • Ej_^ where a is the polarizability. 
In Figure 9 it is shown how a polarizable medium screens out the electric field 
and decreases the work function. These considerations apply for dielectrica 
and metals, if they are not in contact with the substrate. The strong vertical 
distance dependence of the electric field in the surface dipole layer involves 
a correlation with the screening induced work function shift A^^. In the 
case of contact of a metallic overlayer with the Fermi level of the substrate, 
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Figure 9: Schematic view of the effect of a thin insulated, polarizable single 
layer (SL) on the work function. The layer gets polarized and accordingly 
reduces the surface dipole. The dashed line is the electrostatic potential 
without single layer, while the solid line shows the effect of the single layer. 



no vertical dipoles in the sense of Figure 9 are induced, although the work 
function may change. However, as we know from the Smoluchowski effect 
[54] a corrugation induces a non-uniform surface charge density and thus 
lateral electric fields (see next Section). 

For a dielectric, as it is h-BN, the screening causes, e.g. in the case of 
h-BN/Ni(lll) a work function lowering of 1.6 eV [36], which corresponds to 
an induced dipole of 0.3 Debye. In this case charge gets displaced, but not 
transferred from the sp^ layer to the substrate. If the medium is metallic, 
as it is for most graphene cases, we expect a charge transfer that aligns the 
chemical potentials of the substrate and the overlayer. A strong interaction 
also alters the chemical potentials and in turn influences the charge transfer. 
In any case the charge redistribution changes the work function, i.e. the 
surface dipole. For weakly bound graphene theory predicts a correlation 
between the work function of the substrate and the doping level [55] . 

2.3.2 Corrugated layers: Lateral electric fields, dipole rings 

If the surface is flat the lateral electric fields will only vary for ionic compo- 
nents in the two sublattices of the sp^ networks. For the case of h-BN the 
different electronegativities of the two elements cause a local charge transfer 
from the boron to the nitrogen atoms. By means of density functional theory 
calculations it has been found that for a free-standing h-BN sheet about 0.56 
electrons are transferred from B to N (B,N)=(0.56,-0.56)e~ [36]. This yields 
a sizable Madelung contribution to the lattice stability of 2.4 eV per 1x1 unit 
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Figure 10: Contour plot of the electrostatic potential at the surface of a 
corrugated single layer on a flat metal. The energies continuously increase 
from the Fermi level towards the vacuum level Ey at z = oc, with 
Ef < El < E2 < Es < E4 < Ey. Left for a single layer dielectric, right for 
a single layer metal. 

cell (see Appendix A. 2). If the h-BN sits on a metal this Madelung energy 
is reduced by a factor of 1/2 due to the screening of the ionic charges. For 
h-BN on Ni(lll) the ionicity slightly increases (B,N) = (0.65,-0.59) e~, where 
the net charge displacement of 0.06 e~ towards the substrate is consistent 
with the above mentioned work function decrease due to polarization [36] . It 
is expected that the relatively strong local ionicity of the h-BN surface has 
an influence on the diflFusion of atoms with the size of the BN bond length of 
0.15 nm. If the surface is not flat, but corrugated, this causes lateral electric 
fields on the length scale of the corrugation, which is particularly interesting 
for trapping atoms or molecules [17]. Lateral fields may occur in dielectric 
or above metallic overlayers. Figure 10 shows the lateral and vertical elec- 
trostatic potential in a corrugated single layer nanostructure. For the case of 
a dielectric (Figure 10 a)) a different distance of the layer from the surface 
imposes a different screening and lateral potential variations, also within the 
dielectric. For the case of a metal (Figure 10 b)) a corrugation causes lat- 
eral potential variations, reminiscent to the Smoluchowski effect [54], but not 
within the layer. Also, it was found for g/Ru(0001) that the metallic case 
causes smaller lateral potential variations [50] . 

For dielectrics or insulators the key for the understanding of the lateral 
electrostatic potential variation came from the a band splitting (see Section 
2.4.2). This splitting, which is in the order of 1 eV, is also refiected in Nls 
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core level x-ray photoelectron spectra [56], where the peak assignment is in 
line with the a band assignment [49, 48]. Without influence of the substrate 
the energy of the a band, which reflects the in plane sp^ bonds, are referred 
to the vacuum level. This means that the sum of the work function and the 
binding energy as referred to the Fermi level, is a constant. Vacuum level 
alignment arises for physisorbed systems, as e.g. for noble gases [57, 58], 
or as it was proposed for h-BN fllms on transition metals [59]. The a band 
splitting causes the conceptual problem of aligning the vacuum level and the 
Fermi level with two different work functions. The local work function, or 
the electrostatic potential near the surface may, however, be different. This 
electrostatic potential with respect to the Fermi level can be measured with 
photoemission of adsorbed Xe (PAX) [60, 57]. Xenon does not bond strongly 
to the substrate and thus the core level binding energies are a measure for 
the potential energy difference between the site of the Xe core (Xe has a van 
der Waals radius of about 0.2 nm) and the Fermi level. Recently, the method 
of photoemission from adsorbed Xe was successfully applied to explore the 
electrostatic energy landscape of h-BN/Rh(lll) [17]. In accordance with 
density functional theory calculations it was found that the electrostatic po- 
tential at the Xe cores in the holes of the h-BN/Rh(lll) nanomesh is 0.3 
eV lower than that on the wires. This has implications for the functionality, 
since these sizable potential gradients polarize molecules, and it may be used 
as an electrostatic trap for molecules or negative ions. The peculiar elec- 
trostatic landscape has been rationalized with dipole rings, where in plane 
dipoles, sitting on the rim of the corrugations produce the electrostatic po- 
tential well. For in plane dipoles that sit on a ring the electrostatic potential 
energy in the center of the ring becomes 

where e is the elementary charge, R the radius of the hole and P = |pj 
the sum of the absolute values of the dipoles on the ring. For R—1 nm and 
A£'pot=0.3 eV, P gets 10 Debye, which is equivalent to 5.4 water molecules 
with the hydrogen atoms pointing to the center of the holes. The strong 
lateral electric flelds in the BN nanomesh may be exploited for trapping 
molecules, or negatively charged particles. It can also act as an array of 
electrostatic nanolenses for slow charged particles that approach or leave the 
surface. Figure 11 shows the electrostatic potential in a dipole ring and 
the square of the related electric fields. The square of the gradient of the 
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Figure 11: Schematic drawing of a dipole ring where in plane dipoles are 
sitting on a ring with radius R. a) the electrostatic potential energy Epot and 
b) the polarization energy, which is proportional to {VEpotY- From [43]. 

electrostatic potential, or the electric field, is expected to be proportional to 
the polarization induced bond energy Epof. 

Epoi = E|| . p,nd = a-El (4) 

where Ey is the lateral electric field, Pind the induced dipole, and a the polar- 
izability. For the case of h-BN/Rh(lll) the origin of the in plane electrostatic 
fields is not the Smoluchowski effect, where the delocalisation of the electrons 
at steps forms in plane dipoles [54]. It is due to the contact of differently 
bonded boron nitride with different local work functions since the screening 
depends on the vertical displacement of the dielectric from the metal, i.e. on 
the corrugation (see Section 2.3.1). 

The concept of dipole rings is not restricted to the above mentioned po- 
larization of corrugated dielectrics. As it was recently shown for the case of 
g/Ru(0001), lateral electric fields also occur above corrugated metals, where 
the dipoles are created due to lateral polarization like in the Smoluchowski 
effect [50]. 

2.4 Electronic structure II: Bandstructure, Fermisur- 
faces 

The band structure of sp^ layers on transition metals has the same signature 
as free-standing single layers. The binding energy difference between the in 
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plane g bands and the out of plane tt bands does not remain constant, when 
the layers come into contact with a substrate. This is a consequence of the 
bonding via the tt orbitals, while for the a bands vacuum level alignment is 
observed. 

2.4.1 Flat layers: The generic case 

The electronic bandstructure of h-BN on Ni(lll) is similar to that of a free- 
standing layer of h-BN, with g and tt bands that have the generic structure of 
a s'p^ honeycomb lattice (see Appendix). In Figure 12 angle resolved photoe- 
mission data from h-BN/Ni(lll) along the VK azimuth are compared with 
density functional theory calculations. In the bottom panel the bandstruc- 
ture calculations are shown for h-BN/Ni(lll) and free-standing h-BN. The 
band width of the g bands increases less than 1 % by the adsorption of the 
h-BN. The widths of the spin split tt bands increase in the spin average by 
4%, where the width is larger for the minority spins. This band structure 
with G and tt bands is generic for all s'p^ systems. Though, the g and the tt 
band do not shift equally in energy upon adsorption of the sy^ layer. The 
shift between the g and the tt bands of about 1.3 eV indicates that the in- 
teraction between the Ni substrate and the h-BN tt and g bands is not the 
same. Also other substrates like Pt(lll) and Pd(lll) cause a similar h-BN 
band structure, although h-BN/Pt(lll) and h-BN/Pd(lll) are systems with 
large negative lattice mismatch [21, 44]. For these 'flat' cases the bonding to 
the substrates is rather independent of the positions within the unit cell. In 
the Section on corrugated layers (2.4.2) we will see that this changes when 
substrates with mismatch and stronger lock-in energy are used. For Ni(lll), 
Pd(lll) and Pt(lll) it was found empirically that the g bands have, within 
± 100 meV, the same binding energy, if they are referred to the vacuum 
level [59] . This vacuum alignment confirms that the g bands do not strongly 
interact with the underlying metal. The binding energy of the tt band, on 
the other hand, varies, and indicates the tt bonding to the substrate. 

For graphene the electronic structure is also well understood [21]. The tt 
bands are of particular interest since they are decisive for the issue on whether 
the overlayer is metallic or not. On substrates the ideal case of free-standing 
graphene may generally not be realized, since the site selective "p^ interaction 
with the substrate is causing a symmetry breaking between the two carbon 
sublattices Ca and Cb- This opens a gap at the Dirac point. The Dirac 
point coincides with the K point of the Brillouin zone of graphene, where 
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the bandstructure is described by cones on which electrons behave quasi 
relativistic (see Appendix). 

The magnitude of the gap is a measure for the anisotropy of the interac- 
tion of the Ca and the Cb carbon atoms with the substrate. The position 
of the Fermi level with respect to the Dirac energy is decisive on whether 
we deal with p-type or n-type graphene (or h-BN). Since the center of the 
gap of h-BN lies below the Fermi level, this means that e.g. h-BN/Ni(lll) 
is n-type. Figure 13 shows the tt band electronic structure of sp^ honey- 
comb lattices near the Fermi energy and around the K point. If the Fermi 
energy lies above the Dirac energy Eoirac-) which is the energy of the center 
of the gap, the majority of the charge carriers move like electrons (n-type). 
Correspondingly, if the Fermi energy lies below the Dirac point the charge 
carriers move like holes (f>-type). This n-type, p-type picture is analogous to 
the semiconductors, where the Fermi level takes, depending on its position, 
the role of the acceptors (p-type) and the donors (n-type) , respectively. Na- 
gashima et al. first measured the tt gaps at K for the case of one and two 
layers of graphene on TaC(lll), where they found n-type behavior and a gap 
of 1.3 eV for the single layer and about 0.3±0.1 eV for the double layer [61]. 

2.4.2 Corrugated layers: a band splitting 

The super structures as described in Section 2.2.2 are also reflected in the 
electronic structure. Figure 14 shows valence band photoemission spectra for 
three different /i-BN single layer systems in normal emission. We would like 
to draw the attention on the a bands at about 5 eV binding energy (see Table 
1). Figure 14 shows valence band photoemission data for h-BN/Ni(lll), h- 
BN/Pd(lll) and h-BN/Rh(lll). Hexagonal boron nitride has 3 a bands 
and one tt band that are occupied. Along F, i.e. perpendicular to the sp^ 
plane, the two low binding energy a bands are degenerate, while the third is 
not accessible for He 1^ radiation. Accordingly the a bands of h-BN Ni(lll) 
and of h-BN/Pd(lll) show one single peak. Those of h-BN/Rh(lll) are 
weak and split in a (7^ and a cr^ contribution [49] . In measuring the spectra 
away from the F point it is seen that both, the cTq. and the bands split in 
two components each, as it is known from angular resolved measurements of 
h-BN single layers with bands only [59] (see Figure 12). 

The (J band splitting indicates two electronically different regions within 
the h-BN/Rh(lll) unit cell. They are related to strongly bound and weakly 
bound h-BN. Later on, h-BN/Ru(0001) [49] was found to be very similar to 
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Figure 12: Experimental gray scale dispersion plot on the TK azimuth from 
angle resolved He la photoemission. The gray scale reproduces the intensity 
(white maximum). Bottom: Spin resolved theoretical band structure of h- 
BN/Ni(lll) along MT and TK in the hexagonal reciprocal unit cell (Brillouin 
zone), b) Spin up, the radius of the circles is proportional to the partial 
Px = Py charge density centered on the nitrogen sites (a bands). Thick yellow 
lines show the bands of a buckled free-standing h-BN monolayer, which have 
been aligned the /i-BN/Ni(lll) a band at F. c) Spin down, the radius of the 
circles is proportional to the partial Pz charge density on the nitrogen sites 
(tt bands). From [36]. 
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Figure 13: Schematic diagram showing the tt band electronic structure of 
sp^ honeycomb lattices around the K point. For free-standing graphene the 
TT and the tt* band lie on cones intersecting at the Dirac point. Interaction 
with a support causes a shift of the Fermi energy, and, if it is anisotropic a 
gap opens. Here {EF^rmi > EDirac)) the case for n-type conduction is shown 
(for details see text). 
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Figure 14: He 1^ normal emission photoemission spectra of h-BN/Ni(lll), 
h-BN/Pd(lll) and h-BN/Rh(lll). While h-BN/Ni(lll) and h-BN/Pd(lll) 
show no sizable a band splitting, the splitting into (7^ and for h- 
BN/Rh(lll) is about leV (From [43]). 
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Table 1: Experimental values (in eV) for the photoemission binding energies 
referred to the Fermi level of the a bands for h-BN single layers on three 
different substrates along F, and the corresponding work functions $. The 
binding energies with respect to the vacuum level are determined by 
E^ + ^. 







Substrate 




E^B'yp 




eW 


Reference 


Ni(lll) 
Rh(lll) 
Pd(lll) 


5.3 
4.57 
4.61 


5.70 


3.5 
4.15 
4.26 


8.8 
8.7 
8.9 


[62] 
[49] 
[43] 



h-BN/Rh(lll). Theoretical efforts [47] and atomically resolved low temper- 
ature STM [48] showed that the h-BN/Rh(lll) nanomesh is a corrugated 
single sheet of h-BN on Rh(lll) (see Section 2.2.2). The peculiar structure 
arises from the site dependence of the interaction with the substrate atoms 
that causes the corrugation of the h-BN sheet, with a height difference be- 
tween strongly bound regions and weakly bound regions of about 0.05nm. 

Interestingly, no a band splitting could be found for the g/Ru(0001) sys- 
tem, although these graphene layers are also strongly corrugated [50]. This is 
likely related to the inverted topography (see Figure 8) and the metallicity of 
graphene on ruthenium. The metallicity can be investigated locally by means 
of scanning tunneling spectroscopy. It has however to be said that tunneling 
spectroscopy is most sensitive to the density of states at the F point but the 
graphene conduction electrons reside around the K point. This experimen- 
tal limitation does not exist in angular resolved photoemission, where on the 
other hand the sub-nanometer resolution of scanning probes is lost. With 
photoemission it is possible to determine the average Fermi surface of the 
sp^ layers. The Fermi surface indicates the doping level, i.e. the number of 
electrons in the conduction band. In the case of graphene cuts across Dirac 
cones (see Figure 13) are measured. If the Fermi surface map is compared 
to another constant energy surface near the Fermi energy, it can be decided 
whether the Dirac point lies below or above the Fermi energy, i.e. whether 
we deal with n- or p-type graphene. In Figure 15 the Fermi surface map as 
measured for g/Ru(0001) and h-BN/Ru(0001) are compared. For the case of 
g/Ru(0001) extra intensity at the K points indicates Dirac cones. The cross 
sections of which correspond to 5% of the area of the Brillouin zone, or 0.1 
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e~. The cross section shown in Figure 15 c) does not resolve a single cone 
with a bimodal distribution around as it is e.g. the case for g/SiC [63], 
but a single peak from which an average cone diameter at the Fermi energy is 
determined. The measurement of the bandstructure along TK indicates that 
these 0.1 e~ are donated from the substrate to the graphene. For the case of 
h-BN/Ru(0001) the Fermi surface shows no peak at K which is expected for 
an insulating or semiconducting sp^ layer [43] . All features that can be seen 
in this Fermi surface map stem from the underlying Ru(OOOl) substrate. 

2.5 Sticking and intercalation 

Graphite is well known for its ability to intercalate atoms [64], and there 
are also reports on intercalation in h-BN [65]. Intercalation is the reversible 
inclusion of guest atoms or molecules between other host molecules. For the 
host material graphite or h-BN intercalation occurs between the honeycomb 
sheets, where the bonding in the sheet is strong and the bonding between the 
sheets is relatively weak. Here, we deal with single sp^ layers on transition 
metals and use the term "intercalation" also for irreversible intercalation, 
as it is observed if metal atoms slip below the sp'^ layer. For the case of 
g/Ni(lll) it has e.g. been found that Cu, Ag and Au intercalate irreversibly, 
although they form no graphite intercalation compounds [30]. 

Intercalation is preceded by sticking (adsorption), and diffusion of the in- 
tercalating species. In the following the model case of cobalt on h-BN/Ni(lll) 
is discussed in more detail. Figure 16 shows the growth of cobalt on h- 
BN/Ni(lll). On flat terraces, as shown in Figure 16 a) three different pat- 
terns are observed, (i) three dimensional (3D) clusters, whose heights scale 
with the lateral diameter, (ii) triangular, two dimensional (2D) islands with a 
constant apparent height and (iii) line patterns [66] . Often these line patterns 
are found to be connected to 2D islands, and 3D islands tend to nucleate on 
such lines. A careful analysis relates the lines with domain boundaries, where 
(B,N)=(fcc,top) and (B,N)=(hcp,top) domains touch [26] (see Section 2.2.1). 
The 2D islands are irreversibly intercalated Co below the h-BN layer, while 
the 3D clusters remain on top of the h-BN. This results e.g. in the prop- 
erty that the 3D islands may be removed, cluster by cluster, with a STM 
manipulation procedure [66]. 

The sticking coefficient, i.e. the probability that an impinging atom sticks 
on the surface, is determined by X-ray photoelectron spectroscopy (XPS). 
While the sample was exposed to a constant flux of about 3.5 monolayers 
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Figure 15: Fermi surface maps, (a) g/Ru(0001). (b) /i-BN/Ru(0001). The 
hexagons indicate the surface Brillouin zones of Ru(OOOl) (red dashed), 
graphite (green sohd) and /i-BN (blue sohd). (c) and (d) show the nor- 
mahzed intensities of azimuthal cuts along the dashed yellow sectors in (a) 
and (b) respectively. From [50]. 
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Figure 16: Cobalt on h-BN/Ni(lll). a) Scanning tunneling microscopy im- 
age, showing triangular two dimensional (2D) intercalated Co islands, and 
circular three dimensional (3D) clusters and a defect line connecting to a 
intercalated island (27x27 nm). b) Schematic side view of the situation, 
showing the topology of the 3D and the 2D agglomerates. Data from [66]. 

of Co per hour, the measurement of the Co uptake on the sample gives a 
measure for the sticking probability [66]. Figure 17 a) shows the temperature 
dependence of the sticking coefficient of cobalt atoms that are evaporated 
by sublimation onto h-BN/Ni(lll). Clearly, the sticking is not unity, as it 
is commonly assumed for a metal that is evaporated onto a metal. Also, 
the sticking is strongly temperature dependent. This means that at higher 
temperature more Co atoms scatter back into the vacuum and that the bond 
energy of the individual Co atoms must be fairly small. The solid line in 
Figure 17 a) indicates the result of an extended Kisliuk model [67] for the 
sticking of Co on h-BN/Ni(lll). This model predicts about 30% of the Co 
atoms not to thermalize on the surface and to directly scatter back, and says 
that the activation energy for diffusion is about 190 meV smaller than the 
desorption energy [66]. 

The intercalation is also thermally activated. In Figure 17 b) the inter- 
calated amount of Co as compared to the total amount of Co on the surface 
is shown for low Co coverages. The fact that at low substrate temperatures 
almost no Co slips below the /i-BN, indicates that the intercalation must 
be thermally activated. The number A^^^ of intercalated atoms, may be 
modeled with a simple ansatz like: 

N^^ = hexp-^ (5) 
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Figure 17: Temperature dependence of a) sticking coefficient of Co from a 
sublimation source at about 1400 K and b) intercalation of Co below h-BN 
on Ni(lll). Data from [66]. 

where k2 is the reaction constant that comprises the sticking and the 
diffusion of Co to the intercalation site, E'j^ the thermal activation energy 
for intercalation at this site and fc^T the thermal energy. For the 3D clusters 
the number N^^ of atoms in the clusters is assumed to be constant: 



where ks is the reaction constant that comprises the sticking and the diffusion 
of Co to the 3D cluster nucleation sites. The backreaction, i.e. the dissolution 
of atoms from 2D islands or 3D clusters is neglected, since there the Co bond 
energies are much larger. The solid line in Figure 17 b) shows a fit with 
E'^ = 0.24 ± O.leV and a ratio |^ = ^fo* -^^ temperature dependence 

of the sticking and the diffusion cancel. A temperature dependence of the 3D 
nucleation site density is not included in this kinetic modeling, because the 
data do not allow the extraction of more than two independent parameters. 
The small value of is an indication that a Co atom finds an intercalation 
site more easily than a 3D cluster. This is consistent with the assignment 
that the defect lines (see Figure 16 a)) act as collectors for the intercalation, 
since the probability for a diffusing particle to hit a line is much larger than 
hitting a point like 3D cluster. 

2.6 Functionality 

The functionality of sp^ single layers bases on a set of attractive properties 
like a, compared to clean transition metal surfaces, low reactivity and high 
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Figure 18: Cqq on h-BN/Ni(lll). a) Scanning tunneling microscopy shows 
that a monolayer of Cqq wets the substrate and forms a regular 4x4 structure, 
b) The valenceband photoemission spectrum is dominated by molecular or- 
bitals of Cqq. The intensity and energy position of the orbitals is strongly 
temperature dependent. Note the shift of the photoemission leading edge of 
80 meV, which corresponds to the transfer of about half an electron onto the 
molecules in going from 150 K to room temperature. From [69]. 

thermal stability. They also feature a change in electron transport perpen- 
dicular and parallel to the surface. In a metal sp^-layer metal hetero-j unction 
the layer changes the phase matching of the electrons in the two electrodes. 
This is particularly interesting in view of spintronic applications where in a 
magnetic hetero-j unction the resistance for the two spin components is not 
the same, and magnetoresistance or spin filtering is expected [68]. If the 
layers are insulating or dielectric, as it is the case for h-BN they act as atom- 
ically sharp tunneling junctions. The corrugated sp'^ layers bear, on top of 
the potential of the flat layers, the possibility of using them as templates for 
molecular architectures. 

2.6.1 Flat layers: Tunneling junctions 

An insulating single layer on a metal acts as a tunneling junction between 
the substrate and the adsorbate. It decouples adsorbates from the metal. 
This decoupling also influences the time scale on which the electronic system 
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equilibrates. The electronic equilibration time scale increases and eventually 
becomes comparable to that of molecular motion. This opens the doors for 
physics beyond the Born-Oppenheimer approximation, where it is assumed 
that the electronic system is always in equilibrium with the given molecular 
coordinates. Such an example is the C6o/h-BN/Ni(lll) junction [69]. Figure 
18 a) shows 3 substrate terraces separated by the Ni(lll) terrace height 
of 0.2 nm. The whole surface is wetted by a hexagonally closed packed 
Ceo layer with a 4x4 super structure, with one Ceo on 4 h-BN/Ni(lll) unit 
cells. The electronic structure of the Ceo layer shows a distinct temperature 
dependence that is in line with the onset of molecular libration. The freezing 
of the molecular rotation at low temperature is well known from Ceo in 
the bulk of fuUerene [70] and Ceo at the surface of fuUerene [71]. The Ceo 
molecular orbitals as measured with photoemission shift by about 200 meV 
in going from 100 K to room temperature. The shift is parallel to the work 
function, which indicates a vacuum level alignment of the molecular orbitals 
of Ceo- This shift signals a significant charge transfer from the substrate 
to the adsorbate. As can be seen in the inset of Figure 18 b) it causes 
an upshift of the leading edge in photoemission. This upshift translates 
in an average charge transfer of about 0.5 e~ onto the lowest unoccupied 
molecular orbital (LUMO) of the Ceo cage. The LUMO occupancy therefore 
changes by about a factor of 7 in going from 100 K to room temperature, 
which is a molecular switch function. The normal emission photoemission 
intensity of the highest occupied molecular orbital (HOMO) as a function 
of the substrate temperature shows the phase transition between 100 K and 
room temperature. The fact that the phase transition is visible in the normal 
emission intensity gives a direct hint that the molecular motion that leads 
to the phase transition must be of rocking type (cartwheel) since azimuthal 
rotation of the Ceo molecules would not alter this intensity. The experiments 
also show that Ceo is indeed weakly bound to the substrate, since it desorbs 
at temperatures of ^500 K, which is only 9% higher than the desorption 
temperature from bulk Ceo- If we translate the work function shift in the 
phase transition into a pyroelectric coefficient pi = dPs/dT^ where Ps is the 
polarization we get an extraordinary high value of -2000 /iC/m^K, which is 
larger than that of the best bulk materials. 

The correlation of the molecular orientation with the charge on the LUMO 
can be rationalized with the shape of the LUMO of the Ceo- The LUMO wave 
function is localized on the pentagons of the Ceo cage and at low temperature 
Ceo does not expose the pentagons towards the h-BN/Ni(lll). Therefore the 



32 



charge transfer is triggered by the onset of rocking motion where the LUMO 
orbitals get a larger overlap with the Fermi sea of the nickel metal. In turn 
this will increase electron tunneling to the LUMO. In an adiabatic picture, 
however, the back tunneling probability is as large as the forth tunneling 
probability, and the magnitude of the charge transfer effect could not be 
explained. It was therefore argued that the magnitude of the effect is a hint 
for non adiabatic processes, i.e. that the back tunneling rate gets lower due to 
electron self trapping. This self trapping can only be efficient, if the electron 
tunneling rate is low, i.e. if the electron resides for times that the molecule 
needs to change its coordinates. 

2.6.2 Corrugated layers: Templates 

Templates are structures that are able to host objects in a regular way. 
In nanoscience they play a key role and act as a scaffold or construction 
lot for supramolecular self assembly that allow massive parallel production 
processes on the nanometer scale. Therefore the understanding of the tem- 
plate function is of paramount importance. It requires the exploration of 
the atomic structure that defines the template unit cell geometry, and the 
electronic structure that imposes the bond energy landscape for the host 
atoms or molecules. On surfaces the bond energy landscape has a deep val- 
ley perpendicular to the surface that is periodically corrugated parallel to 
the surface. The perpendicular valley is responsible for the adsorption and 
the parallel corrugation governs surface diffusion. If the surface shall act as a 
template, and e.g. impose a lateral ordering on a length scale larger than the 
(1x1) unit cell, it has to rely on reconstruction and the formation of super 
structures. For the case of the h-BN/Rh(lll) the superstructure has a size 
of 12x12 Rh(lll) unit cells on top of which 13x13 h-BN unit cells coincide 
(see Section 2.2.2). It is a corrugated sp^ layer which displays a particular 
template function for molecular objects with a size that corresponds with the 
nanostructure. The super cell divides into different regions, the wires, where 
the layer is weakly bound to the substrate and the " holes" or " pores" where 
h-BN is tightly bound to the substrate. The holes have a diameter of 2 nm 
and are separated by the lattice constant of 3.2 nm (see Section 2.2.2). 

Figure 19 documents the template function of h-BN/Rh(lll) for three 
different molecules at room temperature. For Ceo with a van der Waals di- 
ameter of about 1 nm it is found that h-BN is wetted, and 12 Ceo molecules 
sit in one Rh (12x12) unit cell. It can also be seen that the ordering is not 
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Figure 19: Room temperature scanning tunneling microscopy (STM) im- 
ages of molecules trapped in li-BN/Rli(lll) nanomesh. a) Ceo: Individual 
molecules are imaged throughout this region, following closely the topogra- 
phy (15x15 nm). The positions in the hole centers are occupied by either zero 
or one Ceo molecule; at two places, large protrusions may represent additional 
corralled molecules (From [45]). b) Naphthalocyanine (Nc) (C48H2eN8): Site- 
selective adsorption (120x120 nm). The inset (19x19 nm) on the top right 
shows an enlargement. The inset on the right is a schematic representation 
of the molecule in the h-BN nanomesh. (From [48]). c) Copper phthalocya- 
nine (Cu-Pc) (C32HieCuN8) molecules trapped at the rim of the holes (18 x 
18 nm). At the given tunneling conditions, the mesh wires map dark, the 
holes map grey and Cu-Pc are imaged as bright objects. The inset shows a 
magnified model of the trapped Cu-Pc molecule. (From [17]). 

absolutely perfect. There are super cells with one Ceo missing and such where 
one extra Ceo is "coralled" in the holes of the h-BN nanomesh. If we take a 
molecule that has the size of the holes of 2 nm, here naphthalocyanine (Nc) 
(C48H2eN8), we observe that the molecules self assemble in the holes. Appar- 
ently, the trapping potential is larger than the molecule-molecule interaction 
that would lead to the formation of Nc islands with touching molecules. If we 
take a molecule with an intermediate size, i.e. copper phthalocyanine (Cu- 
Pc) (C32HieCuN8) with a van der Waals diameter of 1.5 nm, it is seen that 
they also assemble in the nanomesh holes. There is a significant additional 
feature compared to the Nc case, i.e. Cu-Pc does not sit in the center of the 
holes but likes to bind at the rims. This observation gives an important hint 
to the bond energy landscape within this superstructure. 

In Section 2.3.2 it was shown that corrugated sp'^ layer systems have 
relatively large lateral electric fields due to dipole rings. These fields may 
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Figure 20: Temperature dependence of normal emission valence band pho- 
toemission spectra from Xe/li-BN/Rli(lll) [17] and Xe/g/Ru(0001) [50]. a) 
Energy distribution curves extracted for three different temperatures (66, 75, 
and 82 K) for Xe/li-BN/Rli(lll). b) Spectral weight of the Xe on the wires 
as a function of temperature, c) Spectral weight of the Xe in the holes as a 
function of temperature, d) Energy distribution curves extracted for three 
different temperatures (67, 87, and 96 K) for Xe/g/Ru(0001). e) Spectral 
weight of the Xe on the mounds as a function of temperature, f) Spectral 
weight of the Xe in the valleys as a function of temperature. The solid lines 
in b), c), e) and f) are fits obtained from a zero-order desorption model. 
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/i-BN/Rh(lll) 


g/Ru(0001) 


Phase 














^d(meV) 
Ni 


181 
25 


184 
17 


208 
12 


222 
13 


222 
14 


234 
23 



Table 2: Experimentally determined Xe desorption energies and Xe atoms 
per unit cell at full coverages Ni for /i-BN/Rh(lll) [17] and g/Ru(0001) [50]. 
For all fits an attempt frequency z/ of 1.2 x 10^^ Hz has been used. 

polarize molecules and provide an additional bond energy. Figure 11 shows 
that this bonding scales with a- Ep where E\\ is the lateral electric field, and 
a the polarizability of the molecule. E"^ is largest at the rims of the h-BN 
nanomesh holes or the g/Ru(0001) mounds. This feature in the bond energy 
landscape is in line with the observation that molecules like Cu-Pc like to 
sit at the corrugation rims of the sp^ layers. A quantitative measure for 
the adsorption energy can be obtained from thermal desorption spectroscopy 
(TDS) which was performed for xenon. 

Figure 20 shows a comparison of Xe/h-BN/Rh(lll) and Xe/g/Ru(0001) 
TDS, where the surfaces are heated with a constant heating rate (3 = dT/dt^ 
and where the remaining Xe on the surface was monitored with photoemission 
from adsorbed xenon (PAX). Since the Xe core level energy is sensitive to 
the local electrostatic potential it also indicates where it is sitting in the sp^ 
super cell. This allows a TDS experiment where the Xe bond energy on 
different sites in the super cell is inferred [17, 50]. 

The desorption data for Xe/h-BN/Rh(lll) on the wires and in the holes 
are shown in Figure 20 b) and c), respectively, and for Xe/g/Ru(0001) on the 
mounds and the valley in Figure 20 e) and f). The data of the weakly bound 
sp'^ layer regions (wires and mounds) are well described with a zero order 
desorption model where dN molecules desorb on the temperature increase 
dT: 

_«=^.exp(-^).dr (7) 

where u is the attempt frequency in the order of 10^^ Hz, (3 the heating 
rate, E^ the desorption energy and ksT the thermal energy. The data are 
fitted to the integrals of Equation 7 where the initial coverage A^i is taken 
from the intensity of the photoemission peaks, the sizes of the Xe atoms and 
the super cells. From this it is e.g. found that 25 Xe atoms cover the wires 
in the h-BN/Rh(lll) unit ceU. The desorption energies of 169, 181, 222 
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and 249 meV for Xe/Xe [72], Xe^/h-BN/Rh(lll), Xe^/g/Ru(0001) and for 
Xe/graphite [73] indicate that they are similar, but have a trend to increase 
in going to more metaUic substrates. The initial coverages and desorption 
energies are summarized in Table 2. For the strongly bound regions (hole 
and valley) this single desorption energy picture does not hold. For Xe/h- 
BN/Rh(lll) it turned out that 12 Xe atoms in the holes have to be described 
with a larger bond energy. Correspondingly different phases C and R have 
been introduced, where C stands for coexistence and R for ring or rim [17]. 
Twelve Xe atoms fit on the rim of the h-BN/Rh(lll) nanomesh holes and 
thus a dipole ring induced extra bond energy of 13% was inferred. For 
g/Ru(0001) the situation is less clear cut, since a two phase fit fits the data 
less well than those of the h-BN/Rh(lll) case. This may be due to the 
different topography of the two systems (see Figure 8). The fact that the 
extra bond energy in Xe/g/Ru(0001) is 4% only is rationalized with the lower 
local work function difference, i.e. the smaller Xe core level energy splitting 
in Xe/g/Ru(0001) (240 meV) than in Xe/h-BN/Rh(lll) (310 meV). 
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A Atomic and Electronic Structure in Real 
and Reciprocal Space 

Graphene and hexagonal boron nitride layers are isoelectronic. Both, the 
(C,C) and the (B,N) building blocks have 12 electrons and both have a 
honeycomb network structure with fairly strong bonds of 6.3 and 4.0 eV for 
C=C and B=N, respectively. The electronic configurations of the constituent 
atoms are shown in Figure 21. 

A.l sp^ hybridization 

The strong bonding within the sp^ network foots on the hybridization of 
the 2s and the 2p valence orbit als. The 2s and the 2p energies in atomic 
boron, carbon and nitrogen are close in energy (in the order of 10 eV) and 
have similar spacial extension. If the atoms are assembled into molecules, 
where the overlap between adjacent bonding orbit als is maximised, the linear 
combination of the s and p orbitals provides higher overlap than that of two 
2p orbitals. Figure 22 shows a s and a Px wavefunction and their spx hybrid, 
which is a coherent sum of the wavefunction amplitudes. The phase in the 
p wavefunction depends on the direction and produces the lobe along the 
X direction, which allows a large overlap with the wave functions of the 
neighboring atoms. 

The hybrid orbitals are obtained as a linear combination of atomic or- 
bitals. For the case of carbon sp^ ^ sp^ and sp^ hybrids may be formed, where 
CH4 (methene), C2H4 (ethene) and C2H2 (acetylene) are the simplest hydro- 
carbons representing these hybrids. The tetrahedral symmetry of sp^ and the 
planar symmetry of sp^ is also refiected in the two allotropes diamond and 



5 


6 




B 


c 


N 


1s2 2s2 2pi 


1s2 2s2 2p2 


1s2 2s2 2p3 



Figure 21: Boron, carbon and nitrogen are neighbors in the periodic table. 
Accordingly, h-BN and graphene are isoelectronic and are constituted by 12 
electrons per unit cell. 
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|S> + |Px> = |sPx> 

Figure 22: The coherent sum of a 1^) and a \px) wavefunction leads to a 
directional \sp^) hybrid. The phases of the wavefunctions are marked with 
+ and where negative amplitudes are shaded. 

graphite, where graphite is the thermodynamically stable phase. For boron 
nitride, accordingly cubic and hexagonal boron nitrides are found, where here 
the cubic form is thermodynamically more stable at room temperature and 
1 bar pressure. 

The sp hybrid orbitals are obtained in combining s and p orbitals. This 
corresponds to a new base for the assembly of atoms into molecules. In the 
sp^ hybrid orbitals the index 2 indicates that two p orbitals are mixed with 
the s orbital. The 2^ orbital is combined with 2px and 2py orbitals, into 3 
two fold spin degenerate orbitals that form the a bonds: 

(8) 
(9) 
(10) 

These 3 a orbitals contain a mixture of | of an 5 and | of a p electron. It 
can easily be seen that the three orbitals lie in a plane and point in directions 
separated by angles of 120°. \sp\) points into the direction [y^, 0,0], \sp\) 

into the direction [-y^, 0], and from the scalar product between the two 
directions, we get the angle of 120°. In the case of the sp^ hybridization the 
fourth orbital has pure p character and forms tt bonds. 

\spI) = \Pz) (11) 
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\spj^ is perpendicular to the sp^ a bonding plane. For the case of graphene 
these Pz orbitals on the honeycomb sp'^ network are responsible for the spec- 
tacular electronic properties of the conduction electrons in the tt bands be- 
cause they are occupied with one electron. 

A. 2 Electronic band structure 

The sp'^ hybridization determines the atomic structure of both, hexagonal 
boron nitride and graphene sheets. They form a two dimensional honeycomb 
structure as shown in Figure 23. The lattice can be described as superposition 
of two coupled sublattices A and B (see Figure 23 a)). In the case of graphene 
both sublattices are occupied by one carbon atom ((7^, C^), while in the case 
of h-BN one sublattice is occupied by boron atoms and the other sublattice 
by nitrogen atoms (B,N). The interference of the electrons between these 
lattices causes the peculiar electronic structure of sp'^ layer networks. The 
12 electrons in the unit cell are filled into 4 Is core levels, and into the 
16 sp^ hybrids that form 3 a bonding, one tt bonding, one tt* antibonding 
and 3(7* anti bonding bands. In the bonding bands the two adjacent sp'^ 
hybrids are in phase, while in the antibonding case they are not. The atomic 
orbitals sp^^ sp^ and spl constitute the in plane a bands, while the sp\ — Pz 
orbitals form the tt bands. From this it can be seen that for flat layers the 
(7 and the tt electrons do not interfere, i.e. may be treated independently. 
In Figure 23 the real space lattice and the corresponding Brillouin zone in 
reciprocal space are shown. In /c-space the reciprocal lattice constant is ^, 
where a is the lattice constant of graphene or h-BN, of about 0.25 nm. For 
the case of graphene, where the two sublattices are indistinguishable, this 
leads to a gap-less semiconductor with Dirac points at the K points. In h- 
BN, the distinguishability of boron and nitrogen leads to an insulator, where 
the TT valence band is mainly constituted by the nitrogen sublattice, and the 
conduction band by the boron sublattice. 

The different electronegativities of boron and nitrogen lead for free-standing 
h-BN to 0.56 e~ transferred from B to N [36]. This ionicity produces a 
Madelung energy EMad- 

1 

Euad = O^Mad ' * — (12) 

47reo a 

with aMad=l-336 for the honeycomb lattice [74], q the displaced charge (in 
the above case 0.56 e~) and a the lattice constant. It has to be noted that this 
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r-space 



k-space 



Figure 23: Honeycomb structure like that of a single layer h-BN or graphene. 
a) real space: the sp^ hybridization causes the formation of two coupled 
sublattices A and B with lattice vectors ai and a2 and lattice constant a = 
= 1^2 1- b) reciprocal space: Brillouin zone. The high symmetry points 
r, M and K are marked, where the reciprocal distance TK is The 
distinction of K and is possible if the system is three fold symmetric 
(trigonal), but not six-fold symmetric (hexagonal). 

Madelung energy applies for the free-standing case. If the ionic honeycomb 
layer sits on top of a metal, the energy in Equation 12 reduces by a factor of 



A. 2.1 TT bands 

It is instructive to recall the basic statements within the framework of the 
"tight binding" scheme of Wallace that he developed for the band theory of 
graphite [18]. Also the description of tight binding calculations of molecules 
and solids of Saito and Dresselhaus is recommended [75] and the most recent 
review of Castro Neto et al. [76]. Tight binding means that we approximate 
the wave functions as superpositions of atomic Pz wave functions on sub- 
lattice A and 5, respectively. Furthermore, they are Bloch functions with 
the periodicity of the lattice. The essential physics lies in the interference 
between the two sublattices. It is just another beautiful example for quan- 
tum physics with two interfering systems. In order to solve the Schrodinger 
equation with this Ansatz we have to solve the secular equation: 



1/2. 



Haa — E 
Hba 




E 







(13) 
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where the matrix elements Haa and Hbb describe the energy on the sublat- 
tices A and B and most importantly Hab the hybridization energy due to 
interference or hopping of the tt electrons between the two lattices. 

The solutions for the energies E_ (bonding) and (antibonding) are: 

E± = ^ {Haa + Hbb ± {Haa- Hbb? + ^ab?) (14) 

Here we show the simplest result that explains the essential physics. It 
is the degeneracy of the bonding tt and the antibonding tt* band at the K 
point of the Brillouin zone, if the two sublattices are indistinguishable. This 
means that at K the square root term in Equation 14 has to vanish. 

The tight binding ansatz delivers values for Haa-) Hbb and Hab'- 

Haa = Ea (15) 

Hbb = Eb (16) 

where Ea and Eb are the unperturbed energies of the atoms on sublattices 
A an S. The square of the interference term Hab is fc-dependent and gets: 

\Hab\^ = i\b (l + 4cos(/c^a) cos (/c^^ a /a/3) + 4: co^^ {kya / \/?>)^ (17) 

kx and ky are coordinates in /c-space pointing along x and respectively, 
and a is the lattice constant (see Figure 23). describes the hybridiza- 
tion between the sublattices and is proportional to the electron hopping rate 
between two adjacent sites on sublattice A and sublattice B. Basically it 
determines the tt band width, which turns out to be 'i^AB- H hopping within 
the sublattices, e.g. between two adjacent A-sites, is allowed this leads to k- 
dependent corrections in Haa and Hbb and to a symmetry breaking between 
the TT and the tt* band [18]. If the phase of the electron wave functions is 
considered, the hexagonal symmetry is broken and a trigonal symmetry, i.e. 
two distinct K points, K and K' have to be considerd [77, 76]. For equivalent 
sublattices, i.e. Haa — Hbb the tt bandstructure is given by the hopping be- 
tween the two sublattices i.e. by Hab and Hba^ respectively. From Equation 
17 it is seen that Hab vanishes at the K point, i.e. for k = [k^^ ky) = (0, -^)- 
The fact that the electrons in the two sublattices do not interfere if they are 
at the K point of the Brillouin zone, is a direct consequence of the symmetry 
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of the crystal, and does not change if hopping between non nearest neighbors 
is included in the model. It should also be mentioned that these results are 
expected for the electronic structure of any system that forms a honeycomb 
lattice. 

In Figure 24 the generic tight binding band structure of a sp^ hybridized 
lattice with two sublattices is shown for hopping ^ab between the sublattice 
A and B, only. Figure 24 a) presents the case of graphene, i.e. equivalent 
sublattices A and B. The two Pz electrons from the atoms Ca and Cb fill the 
TT band. The highest energy corresponds to the Fermi energy and lies at the 
K point. This means that the Fermi surface of this system consists of points 
at the K points. The peculiarity that the resulting Fermi surface encloses no 
volume leads to the notion that graphene is a gapless semiconductor. The 
band structure in the vicinity of the K point is most interesting. The dis- 
persion E{k) of the electrons is linear, i.e. (ff )^ = vp — const.. The linear 
dispersion resembles to that of massless photons, or relativistic particles with 
E ^ rrioC^, where rrio is the rest mass. It constitutes so called Dirac cones 
at the K point of the Brillouin zone (see Figure 24) a)). From Equation 17 
we get with H at the K point the Fermi velocity vp = ^'^^ , which is for 
a TT bandwidth of 6 eV or hopping rate of ^ about c/300. In brackets a 
seeming contradiction has to be clarified: The second derivative of E{k) at 
the K point is zero. With the relation for the second derivative of ^ = 
the eflFective mass rrf of the electrons and the holes is infinite, i.e. electrons 
and holes at the K point may not be accelerated, as it is the case for pho- 
tons. However, since the Fermi velocity is not zero, electrons have not to be 
accelerated in order to be transported. 

In Figure 24 b) the result for a lattice with two inequivalent sublattices 
that corresponds to the case of h-BN are shown. For the sake of simplicity, 
the same hybridization has been chosen. The symmetry between the A 
and the B lattice is broken, if the energy of the unperturbed atoms Ea and 
Eb are not the same (see Equation 14), which is obviously the case for boron 
and nitrogen. The band structure is similar, but at the K point a gap with 
the magnitude of \Ea — Eb\ opens, and no Dirac physics is expected. 

A. 2. 2 a bands 

The a bands form the strong bonds between the atoms in the x — y plane. 
For flat layers they are orthogonal to the tt bands and can be treated inde- 
pendently. The case is more involved than that of the tt bands since here the 
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r KMr r KMr 

Figure 24: Tight binding tt and tt* band structures of honeycomb lattices as 
shown in Figure 23 along TKMr. The hopping rate Jab/T^ is kept constant 
for both cases, a) Case for two undistinguishable sublattices A and B with 
Ea=Eb. Note the emergence of a gapless semiconductor: If each sublattice 
contributes one electron, the Fermi surface is constituted by the Dirac point 
at K. b) Case for two distinguishable sublattices A and B with Ea > Eb- 
Note: At X a gap with the magnitude Eg = Ea — Eb opens, and with an 
even number of tt electrons the system is an insulator. 

3 atomic orbitals (5, Px and Py) on the two sublattices give rise to 6 bands. 
Also the overlap between the different atomic orbitals gets larger and the s-p 
mixing is a function of the k vector. With the tight binding ansatz similar to 
Equations 14-17 3 a bands were derived [75]. Essentially the secular equation 
is now the determinant of a 6 x 6 matrix leading to 6 bands. At F the lowest 
lying band, ctq, has s character and the two remaining a bands, ai and a2 
are degenerate and mainly Px and Py derived. 

It is interesting to note that also for the a bands the bandstructure forms 
cones, if the base of the honeycomb lattice is homonuclear (graphene). They 
are reminiscent to the Dirac cone, where the tt and the tt* band touch. At 
K and a binding energy of about 13 eV the ctq and the ai band touch. 
For heteronuclear bases in the honeycomb (h-BN) also a gap opens, as it is 
observed for the tt and tt* band. Of course, the conical band touching is less 
important for graphene and h-BN, since both involved a bands remain fully 
occupied. The measurement of this gap, nevertheless, would open a way to 
distinguish the a and the tt bonding to the substrate. 

Figure 25 shows a state of the art Density Functional Theory (DFT) 
bandstructure calculation for a single layer graphene and h-BN, respectively. 
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Figure 25: Density functional theory (DFT) bandstructure calculation along 
TKMT for graphene and h-BN single layers. The eigenvalues for discrete 
/c-points are shown as dots. The solid lines connect the eigenvalues in the 
Brillouin zone where the band assignment from the tight binding ansatz with 
a 2s and 2p basis set holds: Three a bands ctq, (Ji and a2 and the two tt bands 
TT and TT* are shown. The three a* antibonding bands are more difficult to 
resolve since they are also mixing in the DFT calculation with 3s and 3p 
contributions, a) Graphene. Note the Dirac cone at K at the Fermi level 
Ef^ where the tt and the tt* bands touch. For the ctq and the ai band a 
similar conical band touching occurs at about -13 eV. b) h-BN. Note the 
band narrowing and the opening of a gap at the K point for the tt-tt* and 
the ao-ai band cones. Calculations by courtesy of Peter Blaha. 
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These calculations consider a much larger basis set than the 2s and the 3 2p 
orbitals as it is done in the tight binding picture. Accordingly more bands 
are found. The dots in Figure 25 are the energy eigenvalues for given k vec- 
tors. The three a, the tt and the tt* band reproduce the tight binding result 
[75]. At higher energies also 3s and 3p orbitals contribute eigenvalues and 
the identification of the antibonding bands becomes difficult. 
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